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Abstract. In the light of the new data on the various neutron and proton electromagnetic form factors 
taken in recent years, we update the dispersion-theoretical analysis of the nucleon electromagnetic form 
factors from the mid-nineties. The parametrization of the spectral functions includes constraints from 
unitarity, perturbative QCD, and recent measurements of the neutron charge radius. We obtain a good 
description of most modern form factor data, with the exception of the Jefferson Lab data on G P E /G P M in 
the four-momentum transfer range Q 2 — 3. ..6 GeV 2 . For the magnetic radii of the proton and the neutron 
we find r p M — 0.857 fm and r% t = 0.879 fm, which is consistent with the recent determinations using 
continued fraction expansions. 
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1 Introduction 

The electromagnetic form factors of the nucleon encode in- 
formation about the structure of the nucleon over a wide 
range of scales. Depending on the kinematical conditions, 
they can be determined most precisely from elastic elec- 
tron scattering using Roscnbluth-separation or in double 
polarization experiments. Only in the last decade has it 
become technically feasible to perform the double polar- 
ization experiments. Thus the data base has been consid- 
erably enlarged since our last dispersion analysis was per- 
formed in 1996 [T]P]. This holds in particular for the neu- 
tron form factors, which can only be measured indirectly 
using deuterium or 3 He targets. For recent reviews, see 
e.g. Refs. [HQ- At low momentum transfer, one obtains 
information about the various nucleon rms-radii, which 
are not only of interest by themselves but e.g. the proton 
charge radius r p E has also to be known to a good accuracy 
to perform precision tests of QED in Lamb-shift measure- 
ments, or conversely, one can use such experiments to pin 
down r p E . For some time, low-energy electron scattering 
and Lamb shift determinations led to different values for 
r E , but this discrepancy has been resolved, see Refs. 5 6 
and references therein. One can also extract vector meson- 
nuclcon coupling constants and eventually study the tran- 
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sition from the non-perturbative regime of QCD to the 
perturbative one. Since the data span the range of four- 
momentum transfers from Q 2 ~ to Q 2 ~ 31 GeV 2 , the 
only model-independent method to analyze these data is 
dispersion theory. The most recent dispersion-theoretical 
analysis of the nucleon electromagnetic form factors dates 
back to almost a decade ago, including fits to space-like 
PP as well as time-like data [2]. For some recent vector- 
meson pole-model approaches to the nucleon electromag- 
netic form factors, see Refs. [7||S]. In view of the new data 
and a new data base collected and discussed in Ref. HJ, 
it seems timely to update the work of Refs. PI2] using 
the data basis of [HJ. It is also interesting to investigate 
whether such a general scheme can lead to the pronounced 
structure around Q 2 — 0.2 GeV 2 found for all four form 
factors in the phenomenological analysis of Ref. [Hj ■ 



2 Formalism 

The electromagnetic structure of the nucleon is parame- 
terized by the Dirac (Fi) and Pauli (F2) form factors for 
the proton and the neutron 

F p/n (Q 2 )=F l s (Q 2 )±Fy(Q 2 ), i = l,2, (1) 

with Q 2 the four-momentum transfer of the virtual pho- 
ton (the photon virtuality). Our conventions are such that 
Q 2 > for space-like momentum transfer. We have ex- 
pressed the nucleon form factors in the isospin basis (S — 
isoscalar, V — isovector) which is most appropriate for 
the dispersive analysis. The experimental data are usually 
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given for the Sachs form factors, which are linear combi- 
nations of Fi and F2 : 
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The analysis of the nucleon electromagnetic form fac- 
tors proceeds most directly through the spectral represen- 
tation given by 1 



1 1"°° 



fi 2 + Q 2 



1,2, 1 = S,V, 



(4) 

in terms of the real spectral functions erf (fi 2 ) = Im F/ (p 2 ), 
and the corresponding thresholds are given by /ig = 3M n , 
piQ = 2M V . The spectral functions encode the pertinent 
physics of the nucleon form factors. In the isovector chan- 
nel, the spectral function is build up by the two-pion con- 
tinuum (including the p-resonance) as given by unitarity 
[TH] plus a series of poles, whose masses and residues are fit 
parameters. In the isoscalar channel, we only have poles, 
where the lowest two are given by the ui and the <f> mesons, 
respectively. We restrict the number of poles in both chan- 
nels by the stability condition of Ref. JT] which requires to 
use the minimum number of vector meson poles necessary 
to fit the data. Note that this condition was also used in 
the earlier fits [ECd - Enforcing the correct normaliza- 
tion conditions for all form factors, the experimental value 
of the neutron charge radius and the superconvergence re- 
lations from the perturbative QCD behaviour of the form 
factors, e.g. Fi(Q 2 ) ~ l/(Q 2 ) i+1 as Q 2 tends to infinity, 
reduces the number of free fit parameters considerably. 
For more details on the form of the spectral functions, the 
stability condition, and the way the various constraints 
are included, see Ref. p. 

It is possible to take all pole masses from physical par- 
ticles (except for one which is determined by the con- 
straints discussed above). As a consequence, we have 2 
(3) free parameters in the isovector (isoscalar) channel if 
we restrict the number of poles to 3 (4). Furthermore, we 
have one more free parameter that characterizes the onset 
of the leading logarithms from perturbative QCD. A dif- 
ference to the earlier fits is that the constraint from 
the neutron charge radius has somewhat changed. While 
in an electron-neutron scattering length of b ne = 

(—1.308 ±0.05) • 10 -3 fm was used, the reevaluation of the 
data lead to b ne = (-1.33±0.027±0.03) • KT 3 fm for scat- 
tering off 208 Pband b ne = (-1.44±0.033±0.06)-10- 3 fm 
for 209 Bi ^21- We have performed fits with both values 
but only show results for the Pb value which seems to be 
favored by the G E data. 



1 Note that we work with unsubtracted dispersion relations. 
Since the normalizations of the various form factors are known, 
one could also work with once-subtracted dispersion relations. 



3 Results 

Before discussing the results, we must specify the data to 
which we fit. We use the data basis collected and specified 
by Friedrich and Walcher in Ref. [5]. It consists of a total of 
190 data points, about a quarter of which were not avail- 
able in 1996 (mainly for the neutron). Furthermore, the 
data basis has also been pruned for inconsistent data and 
therefore does not make use of some of the data that were 
included in the fits of For a detailed list of the data 
taken into consideration, we refer the reader to Ref. j^]. 
In particular, the Jefferson Lab data for G P E /G P M [TMH] 
are treated as data for G P E in |5] since G V M is supposed to 
be known better for these virtualities. 

Fitting these data with the readjusted constraints leads 
to the following parameters. The isoscalar masses are M u = 
0.782 GeV, M = 1.019 GeV, M s , = 1.65 GeV, and 
Mg/i = 1.68 GeV. Note that in contrast to the earlier 
fits, we can not work with 3 isoscalar poles only. The 
isovector masses are M p > = 1.05 GeV, M p » = 1.465 GeV, 
and Mpin = 1.70 GeV. Note that except from M p >, whose 
value is fixed from the various constraints, all these masses 
correspond to physical particles listed in the PDG ta- 
bles. The corresponding residua for the isoscalar poles are: 
= 0.767, 0% = 0.318, af = -0.832, af = -0.250, 



2.09, af = 3.97, a 



s" 



2.0-1 . and 



-3.76, 



where the subscript 1 (2) refers to the vector (tensor) 
coupling of the corresponding vector meson to the nu- 
cleon. Similarly, we have for the three isovector poles: 
a( = -0.154, a( = -0.306, af = 1.15, af = -4.42, 

a p = —1.32, and a% = 3.18. The QCD parameters (for 
definitions, see [T]) are 7 = 2.148, A 2 = 18.28 GeV 2 and 
Qo = 0.35 GeV 2 . The value of A 2 characterizes the on- 
set of the leading logarithmic behavior from perturbative 
QCD. We note that the value for A 2 chosen here is about 
a factor of two larger than in the earlier fits PGj- This 
change is due to the strong deviation of the new Jefferson 
Lab data for G E /G P M |13U14j from the asymptotic predic- 
tion of perturbative QCD. 

The errors in the fit parameters quoted above are un- 
known. Due to the highly nonlinear nature of the problem 
(remember that we fit all 4 form factors simultaneously) 
a trustworthy error analysis is a very non-trivial problem. 
This is complicated by the fact that the errors of some 
inputs (e.g. the two-pion continuum |10j) are not known. 
For this reason, no errors were given in the previous dis- 
persion analyses by Holder et al. ^H] and Mergell et al. 
[I] (and also not in the widely used dispersion theoretical 
analysis of pion-nucleon scattering by Hohler and cowork- 
ers QUI)- At present, it is an open problem how to assign 
a serious theoretical error to such type of analysis. 

In the previous dispersion analyses 012] , a strong fine 
tuning of the residues of the excited p mesons in the isovec- 
tor channel could be observed. It is well known that the 
experimental dipole behavior of the form factors requires 
two narrow structures with opposite signs close to each 
other. The fits in the previous analyses, however, tended 
to favor very large residues (~ 40) for a small overall gain 
in x 2 - Taken at face value, these residues implied unphys- 
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Fig. 1. The nucleon form factors for space-like momentum transfer. The solid line gives our best fit for b„ e = ( — 1.33 ± 0.027 ± 
0.03) ■ 10" 3 fm [12 while the dash-dotted line gives fit 1 from Ref. |2|. The dashed lines give the result of the phenomenological 
fit of Ref. 0j. 



ically large couplings for the isovector poles to the nu- 
cleoli. In order to avoid this problem, we have restricted 
the magnitude of the vector meson coupling constants in 
the present analysis. As a consequence, the residues of all 
poles in the isoscalar and isovector channels are of order 



In Fig. 2] we show the resulting Sachs form factors 
by the solid line. The form factors are normalized to the 
canonical dipole fit, with the exception of the neutron elec- 
tric form factor, which is not scaled. We can describe most 
of the data fairly well, with the exception of the Jefferson 
Lab data on G P E in the Q 2 range from 3 to 6 GeV 2 (see 
also the discussion in Refs. |T7ITB] h The x 2 /datum of this 
fit is 2.07. For comparison, we exhibit fit 1 of Ref. [2] by 
the dash-dotted line. This fit is based on the 1996 data 
basis and describes the data clearly worse than our new 
fit. Also shown by the dashed line is the fit with the "phe- 
nomenological ansatz" of |Hj , which has 6 free parameters 



for each of the 4 form factors. 2 It does significantly better 
for G P E . In Ref. [5], a number of other fits with a "phys- 
ically motivated ansatz" were also performed. These fits 
are generally comparable in quality to the one shown in 
Fig. ^ but give a better description of G P M . 

In the fits of Ref. a pronounced structure around 
Q 2 ~ 0.2 GeV 2 was found for G E (and could also be iso- 
lated in the other form factors as a small effect). This re- 
gion of momentum transfers is shown in Fig. [21 in greater 
detail. The structure was attributed to a contribution from 
a very long range pion cloud (for a discussion of this point, 



2 We have reproduced this fit from Eqs. (4-6) and the num- 
bers given in Table 2 of the published version of Ref. For 
G P M , we have used the following values for the parameters: 
oio = 1.0024, an = 0.7554, a 20 = 1.00 00 - o l0 = -0.0024, 
021 = 14.68, a b = -0.159, Q b = 0.325, a b = 0.231 which differ 
slightly from what is given in Table 2 of the published ver- 
sion of Ref. |21 . With the numbers in that table, one can not 
reproduce the curve shown in Fig. 4 of 0. 
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Fig. 2. The neutron charge form factor Gjj for momentum 
transfers Q 2 = 0...0.9 GeV 2 . The curves are as in Fig. Q 



see also Ref. ^Hj). Our analysis does not show the pro- 
nounced structure around Q 2 ~ 0.2 GeV 2 . What this re- 
sult means remains to be seen. Taking the errors of the 
experimental data into account, the structure is not an 
unambiguous feature of the data (cf. Fig. EJ. Ultimately, 
the question of whether the structure found in Ref. 9 is 
real physics or not should be decided by more accurate 
data in the region Q 2 ~ 0.2 GeV 2 . A similar plot for the 
proton charge form factor in the low momentum transfer 
region is shown in Fig. |3J The data in this range are suf- 
ficiently spread to leave room for various interpretations. 



For the electric and magnetic radii of the nucleon, we 
find the values r p E = 0.848 fm, r p M = 0.857 fm, and = 
0.879 fm. As in the previous analysis the proton 

charge radius is somewhat small compared to the recent 
precise determinations r p E — 0.880(15) fm 5 and r p E = 
0.895(18) fm |B] from low-momentum-transfer data and 
r p E = 0.883(14) fm j20j from Lamb shift measurements. 

Note that the recent determinations of the proton charge 
radius by Rosenfelder and Sick are based on analyz- 
ing low-momentum cross section data directly and include 
corrections from two-photon exchange. Our analysis relies 
on form factor data extracted in the one-photon exchange 
approximation where the two-photon corrections enter in 
the systematical error. A comparison of our value to the 
analyses of low-momentum transfer cross sections directly 
allows one to estimate systematic errors in the extraction 
of the radii. Of course, one might contemplate a fit to 
cross section data only. However, for consistency the two- 
photon corrections have to be applied to the full data basis 
(as opposed to low Q 2 data only) and this will certainly 
take a long time. This issue has to be addressed in the fu- 
ture but requires a large common effort by experimenters 
and theorists. 

Our results for the magnetic radii of the proton and 
neutron are consistent with the recent values from con- 
tinued fraction expansions: Kubon et al. |2l extracted 
r M = 0.873(11) fm from precise data for G%[, while Sick 
[Hll22| obtained = 0.855(35) fm from a careful analysis 
of the world data on elastic electron-proton scattering. To 
simplify the comparison, we have collected our results and 
the other recent radius determinations in Tabled 
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Fig. 3. The proton charge form factor G P E for momentum 
transfers Q 2 = 0...1 GeV 2 . The curves are as in Fig. Q 
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Table 1 . Comparison of our results for the radii 



and 



rjjf with other recent determinations. The numbers in paren- 
theses indicate the error in the last digits. 



4 Outlook 



The luNN and <f>NN coupling constants derived from 
the residua given above are gf NN = 21.4, g% NN = 0.9, 
gf NN = —10.3, and g% = —3.1. The absolute values 
of these coupling constants are comparable with the re- 
sults of Refs. but the small tensor couplings g% NN 
and g% NN have changed sign. This change of sign indicates 
that the error in the extracted vector meson coupling con- 
stants is of the order of the tensor couplings. 



We have shown that a dispersion-theoretical analysis based 
on a minimal number of poles can describe most of the cur- 
rent world data on the nucleon electromagnetic form fac- 
tors for space-like momentum transfer. While the charge 
radius of the proton is still somewhat small, the magnetic 
radii of the proton and neutron are in good agreement 
with the recent determinations using continued fraction 
expansions P21II22| . 
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The spectral functions used here can be improved. Ide- 
ally, one would like spectral functions consisting of vari- 
ous continuum contributions with vector mesons emerging 
naturally as resonances and avoid explicit pole terms alto- 
gether. In particular, the isoscalar region around masses of 
1 GeV could be supplemented by explicit irp j22j and K K 
|24| continuum contributions. The <j> meson would then 
appear as a finite width resonance in the KK (and possi- 
bly 7rp) continua, thereby eliminating the need to include 
it explicitly as a pole term. It would also be desirable to 
replace the poles with masses in the region from 1.5 to 2.0 
GeV with an appropriate continuum contribution. How- 
ever, the analytical continuation of experimental data that 
is required to obtain the continua becomes more and more 
difficult the higher one gets in mass, so that the latter goal 
is not realistic in the near future. Another improvement 
is to construct a better representation of the perturbative 
QCD behaviour for large Q 2 . This would allow to include 
all existing time- like data in the fits, and thus lead to a 
more consistent description of these fundamental quanti- 
ties. 

Furthermore, the proton charge radius extracted from 
low-momentum-transfer data and Lamb shift measurements 
could be included as a further constraint in the analy- 
sis. This would allow to check the consistency between 
the high-Q 2 data and the low-Q 2 extractions of the radii. 
Whether this requires the introduction of additional vec- 
tor meson poles is an open question. Last but not least, 
a full error analysis of the extracted radii and coupling 
constants should be carried out. As mentioned above, this 
is a nontrivial task because of the highly nonlinear nature 
of the problem. Work along these lines is under way j2S]. 
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